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■ We present a study of the SU (3) glueball spectrum for all J PC values at lattice spacings 

down to a -1 = 3.73(6) GeV (/? = 6.4) using lattices of size up to 32 4 . We extend previous 
studies and show that the continuum limit has effectively been reached. The number of clearly 
identified J PC states has been substantially increased. There are no clear signals for spin-exotic 
glueballs below 3 GeV. A comparison with current experimental glueball candidates is made. 



Introduction The extraction of reliable predictions for the glueball spectrum of QCD remains 
an important challenge for lattice gauge theory. As part of a recent programme to study non- 
perturbative pure SU(3) gauge theory closer to the continuum limit, we have obtained new data 
for glueball masses which confirm that results of relevance to continuum physics are indeed being 
achieved with currently accessible lattices. The low-lying spectrum for pure glue states below 3 
GeV or so has become well established now and awaits progress with understanding the effects of 
light-quark mixing. 

Glueballs are not particularly light — they start around 1500 MeV — and have no non-trivial flavour 
content. The extraction of a signal in the presence of vacuum fluctuations is therefore more difficult 
than for many other hadrons or for potentials. In this situation it is highly desirable to perform 
coherent measurements over a suitable /3-range, in order not to be lost in possible systematic effects. 
In this spirit we apply here the techniques used by and, in some cases, pioneered by Michael and 
Teper (MT) jL], |2], [| and extend their analysis. They used lattices ranging from 10 4 to 20 4 , at (3 
values up to 6.2. In the meantime, there has been progress both in the available computing power 
and in the efficiency of updating algorithms. In this work we have used a hybrid [Q, || of heat-bath 
and over-relaxation. The code was specifically developed for the Connection Machine and was run 
on an 8K machine at Wuppertal and a 16K machine at Edinburgh. The key parts of the code, 
including the group theory, were thus independent of previous work. We have concentrated on 
P = 6.4 on 32 4 — slightly larger in physical size than the largest size used by MT, but have also 
taken data at (3 = 6.0 and 6.2 where a direct comparison could be made. 

Measurement procedures The (3 = 6.4 results presented in this letter were based on the 
measurement of 3220 configurations, each separated by ten sweeps. Every fifth sweep was a heat- 
bath step, the remainder being Creutz over-relaxation. The data was obtained in two parts, from 
hot and cold starts with at least 2000 sweeps used to equilibrate in each case. During the subsequent 
analysis, described below, a careful check was made that no residual equilibration effects were 
present, that both samples were consistent, and the measurement sampling rate was reasonable 
compared with the autocorrelation times. A direct measurement of the autocorrelation time gave 
r < 20 sweeps for the correlators of interest. In order to allow greater flexibility in analysis and to 
allow further cross-checking, the in-line measurements were done in a rather general way. They were 
made on the (L/2) 3 L lattice configuration obtained from an L 4 configuration by Teper fuzzing || 
with the link/staple mixing parameter a = 1.0. On each time slice, operator momentum transforms 
for a variety of oriented 'shapes' and for all cubic orientations were stored for this level of fuzzing, 
and for each subsequent level up to the maximum physically reasonable. The shapes were as noted 
in table [l] (see ref. [||] for details and a diagram). Non-zero momentum operators were used only 
for the plaquette shape for which we considered k^ = 0, ±1 where = ^-k^. Because of the 
initial fuzzing step before measurement, the space points summed over were spaced by 2 units so 
that the momentum eigenvalues k^ were unique only up to modulo L/2. For small momenta the 
contamination is expected to decay very fast in Euclidean time. By studying similar sized lattices 
in SU(2) we did confirm this effect but found that the t = correlations, and hence 1/0 ratios 



showed significant contamination from the high momentum piece. Because of the variational nature 
of the calculation (see below) this did not affect our spectrum results at all. The total time used 
to update the gauge configuration and make these primary measurements was of order 300 hours 
(16K CM-200 equivalent). The gauge update time was 3.8 //sec per link. 

The operator sums for each time-slice were analysed off-line. We studied Euclidean time correlators 
for all representations of Oh- A\, A2, E, T% and T2 for both values of parity and C-parity |7j. The 
relevant projection table is given, for example, in ref. Q. In addition, we studied Polyakov line 
correlators (the torelon) as a cross check on the string tension. Further data acquired for smeared 
Wilson loops and for the topological susceptibility will be presented elsewhere |J. In the off-line 
analysis of correlations, a variational approach was used Q in which a matrix of correlators is 
formed using, as basis, the different relevant operator shapes and fuzzing levels. By diagonalising 
the transfer matrix and studying ratios of eigenvalues at consecutive Euclidean times, one obtains 
upper bounds on the effective mass (or energy for non-zero momentum) of the ground state in each 
channel. In principle, estimates of excited states can also be made. 

There are two cross-checks on the reliability of the mass values so obtained. First, the overlaps 
for the various operators are obtained. For a stable determination, one would prefer large 'wave- 
function' components carrying the same sign rather than a delicate cancellation (as a result of a 
poor choice of basis). Indeed, we have checked the stability of our results to using smaller and 
differing samples of basis operators. The A~i + receives contributions from a broad range of shapes 
and fuzzing levels, while the remaining states receive dominant contributions from the maximum 
fuzzing level, mostly from the 'hand' shaped loops (see table |]). Second, one expects the ground 
state in each channel to dominate at large Euclidean time. We have monitored the difference 
between successive effective masses to find at which t value this becomes statistically insignificant. 
For the determination of errors, we have always used the bootstrap sampling procedure where the 
data are organised in bins large compared with the measured autocorrelation lengthE For the 
majority of states, the effective mass 'plateau' identified in this way starts at time ratio 2/3. 

Results Table [2] contains the measured effective masses of all ground state glueballs which can be 
studied on a hypercubic lattice. Where a significant signal was found, the chosen plateau value and 
its associated error estimate is indicated by bold face. Where larger time ratios gave higher masses 
or where the plateau was not particularly well established, the error was conservatively estimated 
from the next larger time ratio. In the final column of the table, the glueball masses m are given in 
units of the string tension from Wilson loops. We find that the spectrum proposed by MT from the 
average of their large volume data at /3 = 5.9, 6.0 and 6.2 ( |3], table 7) is in agreement with our 
(3 = 6.4 mass ratios. Our data at 6.0 (on 16 3 x 32) and at 6.2 (on 32 4 ) are themselves consistent 
with the corresponding results of ||. Overall, we conclude that the new lattice measurements of 
the low-lying spectrum at /3 = 6.4 do indeed represent useful information about continuum physics. 

The restoration of symmetry provides an additional and more stringent test for continuum physics. 

1 In fact, for test observables, we have used the bin size dependence of the measured variance to cross-check the 
direct measurements of the autocorrelation time quoted above. 



The expected continuum J content of the Oh representations can be found, for example, in 
ref. (Tj. In table we indicate only the lowest possible J PC . We confirm that the E and T2 ground 
states (that both contribute to J = 2 in the continuum 0(3) symmetry group) exhibit the expected 
degeneracy for all PC combinations. For PC = ++ this has been found previously in ref. 0] for 
(3 > 6.0. A related requirement is the restoration of the continuum dispersion relation i.e. Lorentz 
symmetry. We have been able to test this for the momenta k 2 = 0, 1, 2, 3. A one parameter fit@ of 
E k to 

E k a = sj<* 2 +(^)' (1) 

yields for the data mo« = .425(12) with x 2 /dof = 0.55. The non-zero momentum results for 
the mass, though slightly higher than the zero momentum value given in table §, agree well within 
errors. These two features give strong support to our statement that, for the low-lying states, our 
(3 = 6.4 results are effectively measurements of the continuum glueball spectrum. 

Different continuum J PC states can contribute to a given Oh representation. A priori their level 
ordering is not obvious. In fact we observe that the 1 ++ mass is definitely larger than the 2 ++ 
mass. By assigning the 1 ++ quantum number to the lattice state we have assumed the 'natural' 
ordering m 1 ++ < m 3 ++. However, the data at (3 = 6.4 show the and T±~ + to be approximately 
degenerate. This is not the naive expectation since the lowest contributing J PC values in each 
case are 3 ++ , 6 ++ and 1 ++ , 3 ++ , 4 ++ , respectively. One possibility is that the 3 ++ state is in 
fact lower than the 1 ++ and so gives a common lightest contribution to both lattice states. These 
higher mass states are difficult to observe cleanly at lower (3 where ma is unhelpfully large. Indeed 
this possible degeneracy was not seen by MT at 6.2 or 6.0. Our estimate of m A ++/\/a at 6.4 is 
consistent with the data of ref ||] at 6.2, and with both our 6.2 and 6.0 data. For the + also, 
we have a reasonable signal at 6.4 unlike at 6.2 where an upper limit only was obtained ||. 

According to table § we observe reasonably good signals for 10 states of different continuum J PC 
contents. These are included as the solid circles in fig. |l|. Moreover, we determine upper limits for 
the masses of the remaining 6 states with J < 4. These are also shown in the figure. With our 
lattice resolution (a -1 ~ 3.7 GeV) and statistics we are in a position to trace the signals over larger 
time separations and achieve more stringent upper bounds on masses than previously possible. 
This improves our capability to separate low lying glueball states and establish the spectrum order. 
The 2 ++ is separated by some 6 standard deviations from the lightest (0 ++ ) glueball. Moreover, 
the 2 ^ glueball is found to be significantly heavier than the . Above the 2 , five further 
states have been identified but their ordering cannot yet be determined. 

Before proceeding with further interpretation of our results, it is useful to convert the continuum 
predictions in the last column of table || to an MeV scale. The extra scale displayed in fig. [I] was 
obtained by multiplying the latter numbers by a string tension value of 440 MeV. Direct comparison 
with experiment is valid only within the following two assumptions: (a) that the glueball masses are, 
for some reason, insensitive to light quark mixing where this takes place and (b) that the physical 
string tension estimate (440MeV) 2 based on a model for Regge trajectories is reliable for the pure 

2 For technical reasons, we have not attempted a full correlated error analysis. 




Figure 1: The measured glueball spectrum at f3 = 6.4. Open symbols represent measured upper 
limits. The origin of the MeV scale is described in the text. 



glue sector. Some encouragement for this latter belief is provided by the fact that the scale for light 
hadron masses set in this way is very reasonable. A recent large scale quenched lattice study of the 
light meson and baryon spectrum at similar lattice spacings shows excellent agreement with 
experiment, provided sufficiently small valence quark masses are used for extrapolations. Using 
the p mass to set the scale at (3 = 6.17 yields a -1 = 2.63(4) GeV which is quite consistent with 
2.78(5) GeV 12] and 2.72(3) GeV Q at (3 = 6.2 as obtained from the string tension. Furthermore, 



values of deduced from independent studies of the SU (3) heavy quark potential (256±20 MeV 



[12 1 and 244 ± 8 MeV ||) are compatible with the (unquenched) values found in experiment [13]. 



To set the scale in what follows, we have used string tension values era 2 extracted from Wilson 



loops measured on sufficiently large lattices: 0.168(11) at /3 = 5.7 [15]. 0.073(1) at (3 = 5.9 [JDJ 
0.0476(7) at f3 = 6.0 Q, 0.0251(8) at (3 = 6.2 @ and 0.0138(4) at (3 = 6.4 §. It is interesting to 
note that the latter value is in good agreement with the effective string tension deduced from our 
Polyakov- Wilson line correlator (torelon) 

cr eS a 2 = am tov /L = 0.440(20)/32 = 0.0138(6) . (2) 



One should remember that the latter is subject to a finite size correction of order tt/3L 2 [17] which 
on our size of lattice is small (+0.0010). 

In fig. we show a new compilation of scalar and tensor glueball masses, measured at various 
couplings. Since the lattice corrections to ma/y/aa are expected to be of order a 2 , we display the 
physical masses as a function of a 2 where physical units on both axes have been set as described 
above. The present data is displayed with full symbols and previous data with open symbols 
{(3 = 6.2, 6.0, 5.9 H, p = 5.7 @). We are aware that the string tension results have been obtained 
by slightly different methods but, for present purposes, the ensuing uncertainties are small compared 
to the statistical errors from the glueball masses. It is clear that, at least for the scalar glueball, 
there appears to be little room for uncertainty in any reasonable extrapolation to a 2 = 0. A linear 
fit is shown as an example. The fitted slope to the scalar glueball data suggests a systematic error 
of less than 5% in extrapolation from [3 = 6.4, which would be of the same order as the statistical 
error. Because of this, we henceforth use the results at (3 = 6.4 as an adequate approximation to 
the continuum spectrum. 

Phenomenological considerations With the above scale (y/a = 440 MeV), the scalar glueball 
mass prediction from our /3 = 6.4 data is 

m 0++ = 1550 ±50 MeV (3) 

where the error here is purely statistical. This value is in agreement with previous lattice glueball 
calculations (e.g. refs. @, g, |, |, |J|]). The status of the G(1560) g{| as a 0++ glueball candidate 



has recently been considerably strengthened by its independent observation in pp — > 67 [20] where 
a strong coupling to the 7777, but not to the 7r°7r° channel, is found. The total width is 245 ±50 MeV. 
Clearly, the lattice calculation is quite consistent with this mass. However, the lattice state is not 
far from the broad qq state /o(1400) observed predominantly in tttt but which also couples to r/?] 
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Figure 2: Scalar and tensor glueball masses as a function of a 2 . Full symbols refer to our data. 
Open symbols are taken from refs. [3,17]. Circles are Af + , diamonds E ++ , and squares T 2 ++ . For 
clarity data points at the same (5 values have been separated slightly. 



(and 77). Because of the influence of the /o(975) and the KK channel to which they both couple, 
the width (a few hundred MeV) and indeed the very nature of the /o(1400) is difficult to establish. 
It seems not unlikely that this pure glue state will suffer mixing and be part of a complex system 
involving the above states. Future lattice studies of light quark mixing will be very illuminating 



on this point. Pioneering attempts to study this [21, 22 1 are hampered by the unphysically large 



quark masses currently accessible and the difficulty in acquiring sufficient statistics. 
The above energy scale estimate puts the tensor glueball at 

m 2 ++ = 2270 ± 100 MeV . (4) 

So far, only one experiment has provided evidence of a 2 ++ glueball candidate in this mass range 



[|§. A series of three <fxf> states in the range 2010 to 2340 MeV with widths of 150 to 300 MeV have 
been seen but not yet independently confirmed. The next predicted glueball state, a pseudoscalar 
at around the same mass (table ^), has no suitable experimental candidates currently. The search 
becomes increasingly difficult at high masses where many states overlap and many channels are 
competing. 

As pointed out previously (e.g. ref. ||), the prediction of low- lying exotic states (i.e. non qq quark 
model states) would have interesting theoretical and phenomenological consequences. Michael [ j24[ | 
has recently reviewed the lattice and experimental evidence for glueball and hybrid states with 
these quantum numbers. Our results (table ^) confirm earlier predictions that no exotic glueball 
states are expected below about 3 GeV. On the lattice, each Oh representation corresponding to 
an exotic J : , H , 1 ^ etc. can also receive contributions from higher, but non-exotic, J PC 
so identification is unlikely to be straightforward in the absence of very precise data. The strongest 
exotic signal we have observed is in the Trj 1 channel. There is also some evidence of a signal in the 
channel. These could correspond to a 2^ exotic glueball at around 3.9 ± .7 GeV. However, 
the lowest non-exotic J PC contributing to the would be 3"^ (5 -1 for the ) and so no 
strong conclusion may be drawn. The A\ and T{~ channels also show some sort of signal. These 
do not satisfy the above criteria for plateau identification and so we only quote these as upper 
limits. Experimental confirmation of exotic states in the above mass range is likely to be very 
difficult. 

In conclusion, we have demonstrated that at /3 = 6.4 we are effectively at the continuum limit 
for the quenched glueball spectrum below 3 Gev. To be specific, we observe clear signals for 10 
different continuum states. Thus the ordering of the underlying spectrum is becoming established. 
Further improved studies of lattice glueballs are both practicable and desirable. In the near future, 
machines capable of sustaining 50 to 100 Gflops on QCD will allow a factor of v50 or so reduction 
in statistical errors and hence greatly improved effective mass signals. The present work represents 
a reduction in lattice spacing by 25% and an increase in physical volume by 70% over previous 
studies. Our results show that there is no need to use larger lattices or smaller lattice spacings 
to probe the mass range which is likely to be of most experimental relevance i.e. below 3 GeV. 
However, higher mass states will require larger statistics and lattice spacings such that ma < 1. 
In order to estimate the possible influence of mixing effects due to light quarks it will be vital to 



have increased precision of meson and glueball masses in quenched QCD. This is almost within our 
grasp. 
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Links Shapes Orientations Fuzz levels 
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plaquette 
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1,2,3,4 


6 


rectangle 
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1,2,3 


6 


chair 


12 


1,2,3,4 


8 


hand 


48 


1,2,3 


8 


butterfly 


24 


1,2,3 



Table 1: Glueball operators used. See ref. [3] for a diagram 



O h Rep. 


J PC . . . 


0/1 


1/2 


2/3 


3/4 


4/5 




A++ 


0++ 


0.604(7) 


0.435(8) 


0.415(14) 


0.402(20) 


0.38(3) 


3.52(12) 


A++ 


3++ 


1.552(15) 


1.06(4) 


1.05(13) 


1.0(3) 


— 


8.9(11) 


E++ 


2++ 


0.911(5) 


0.653(11) 


0.620(17) 


0.61(3) 


0.56(6) 


5.25(25) 


rp+ + 


2++ 


0.914(5) 


0.638(9) 


0.598(14) 


0.55(2) 


0.52(4) 


5.07(17) 


T++ 


1++ 


1.657(10) 


1.10(3) 


1.06(8) 


1.0(2) 


0.8(6) 


9.0(7) 


^ + 


o-+ 


1.155(9) 


0.751(18) 


0.63(4) 


0.69(7) 


0.68(14) 


5.3(6) 


A 2 + 


3 -+* 


2.34(3) 


1.56(18) 


2.6(24) 








E~+ 


2-+ 


1.265(8) 


0.853(16) 


0.83(4) 


0.77(8) 


0.9(2) 


7.0(3) 


T 2 + 


2-+ 


1.284(7) 


0.851(12) 


0.79(3) 


0.80(6) 


0.95(19) 


6.7(5) 


Tf+ 


I — h* 


1.824(10) 


1.22(3) 


0.99(11) 


1.2(5) 






Af 


o+-* 


2.24(2) 


1.31(9) 


0.8(2) 


0.9(6) 






Af 


3+- 


2.68(3) 


1.6(2) 


0.9(5) 








E+- 


2 +-* 


2.090(18) 


1.23(5) 


1.2(2) 


0.5(4) 




10(2) 


T+- 


2 +-* 


1.461(12) 


0.97(2) 


0.91(5) 


0.93(12) 


1.0(4) 


7.7(10) 




1+- 


1.188(5) 


0.837(13) 


0.78(3) 


0.82(7) 


0.92(14) 


6.6(6) 


K~ 


— * 


2.24(3) 


1.55(14) 


1.4(6) 










3~ 


2.70(3) 


1.5(2) 










E— 


2~ 


1.715(13) 


1.07(3) 


1.03(9) 


1.7(6) 




8.7(8) 


T 2 - 


2— 


1.804(9) 


1.13(3) 


1.08(9) 


0.9(3) 




9.2(8) 


TT~ 


1~ 


1.845(11) 


1.17(4) 


1.22(13) 


0.8(3) 


1.2(10) 


9.9(11) 



Table 2: Glueball effective masses at (3 = 6.4, in lattice units. The last column shows the ratio to 
the string tension where the quoted errors arise only from the glueball statistical errors. The J PC 
value displayed labels the lowest continuum representation that can contribute. Exotic J PC content 
is indicated by (*). 



